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The dynamical evolution of a quantum system composed of two coupled 
cavities, each containing a two-level atom and a single-mode thermal field, is 
investigated under different conditions. The entanglement between the two 
atoms is controlled by the hopping strength and the detuning between the 
atomic transition and the cavities. We find that when the atomic transition 
is far off-resonant with both the eigenmodes of the coupled cavity system, 
the maximally entangled state for the two atoms can be generated with 
the initial state in which one atom is in the ground state and the other is 
in the excited state. When both the two atoms are initially in the excited 
state, the entanglement exhibits period sudden birth and death. By choosing 
appropriate parameter values, the initial maximal entanglement of the two 
atoms can be frozen. The relation between the concurrence and cooperative 
parameter is calculated. © 2012 Optical Society of America 

OCIS codes: 270.2500, 270.5580 



1. Introduction 

Quantum system inevitably interacts with the surrounding heat environments, i.e., thermal 
fields. The thermal field, which is emitted from a source in thermal equilibrium, has become 
an important topic in quantum theory of damping and quantum experiment of laser [1]. 
Under the condition of constant average energy, the thermal field maximizes the entropy 
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of the field and is always considered to be chaotic [2]. However, Bose et al [3] showed 
that a single two-level atom and a thermal field could be entangled due to their linear 
interaction. Recently, the model in which two atoms interact with a single-mode thermal 
field has attracted great interest [4-10]. Kim et al [4] considered two identical two-level 
atoms which are resonantly coupled to a single-mode thermal field and found the field could 
entangle two atoms which are initially in a separate state, while Zhang [5] generalized Kim's 
idea to the case that two atoms are slightly detuned from the thermal field. El-Orany [6] gave 
a exact treatment for entanglement between two atoms when they interact simultaneously 
with a single-mode thermal field through multiphoton exchange. On the other hand, the 
fundamental dissipative process of the radiation field inside a cavity may be understood from 
a simple model where the mode of interest is coupled to a multi-mode thermal reservoir. There 
have been some studies devoted to the infiuence of thermal noise on the behavior of different 
quantum system [11-16] and to the possibility of the quantum information processing via 
thermal fields [17-19]. However, previous theoretical schemes and experiments focus on the 
case in which two atoms interact with a common single-mode thermal field and the maximal 
entanglement can be obtained only when the atomic transition frequency is highly detuned 
from the field frequency. 

Recently, considerable attentions have been paid to the coupled-cavity system [20-26], 
which is promising to overcome the difficulty of individual addressability existing in a single 
cavity. Dynamics in a coupled-cavity array offers a basic setting for distributed quantum 
information processing and can be studied for a wide range of parameters, such as atom- 
cavity detuning, intercavity coupling rate, etc. In this paper, we investigate the entanglement 
dynamics for two atoms trapped in two coupled cavities, each containing one single-mode 
thermal field. 

Contrary to previous studies in which two atoms simultaneously interact with only one 
single-mode thermal field, here we consider a quantum system comprising two identical two- 
level atoms trapped in two coupled cavities, each of which is initially in a thermal state. The 
entanglement dynamics of the atomic system depends upon the photon hopping strength, the 
detuning between atomic transition frequency and field frequency, and the initial state. We 
find that even when the atoms symmetrically and resonantly interact with the local thermal 
fields, the maximally entangled state for the atoms can be generated with the initial state 
in which one atom is in the ground state and the other is in the excited state. This result is 
obviously different from the case that two atoms symmetrically and resonantly interact with 
a common thermal field [4] , in which the obtainable entanglement is rather low. Under certain 
conditions, the thermal fields can induce the periodical entanglement sudden birth and death 
when the atoms are both initially in excited states or ground states. Remarkably, the peaks 
of the concurrence do not decrease monotonically when the evolution time increases, which 
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is distinguished from previous examples [27-30]. By choosing appropriate photon hopping 
strength and detuning between the atomic transition and the cavities, one can freeze the 
entanglement between the two atoms when they are initially in the maximally entangled 
state, which may have practical application in quantum memory and quantum storage [31- 
33]. The influences of atomic spontaneous emission and photon leakage out of the cavity are 
analyzed. 

2. The Theoretical Model 

The coupled-cavity system under consideration is shown in Fig. 1. Two identical two-level 
atoms 1 and 2 are trapped in two coupled cavities, each containing a single-mode thermal 
fleld p/. = 1, 2). In the frame rotating at the cavity frequency, the total Hamiltonian under 
the rotating- wave approximation is {h = 1): 

H = Ya=i ^ia\ai + gi{Sl ai + o}^ + J {a\a2 + aia\) , (1) 

where = \ei){gi\ and = \gi){ei\ with |e^) and \gi) being the excited state and ground 
state of the ith atom, a] and are the creation and annihilation operators for the ith thermal 
field, gi describes the coupling strength between the ith atom and the ith thermal field, and 
5i is the detuning between the atomic transition and the local cavity mode. J represents the 
coherent photon hopping strength between the two cavities. 

The density operator of the ith thermal radiation field is described by 

Ph = E^=o^^(^^)l^^)(^^L (2) 

= (^ItB^' (3) 

where = (^e^'^Si/^BT^ _ -j^^-i -g ^.j^^ mean photon number and is the Fock state, fc^ is 
the Boltzmann constant, Wf. is the frequency of the cavity mode, and Ti is the temperature. 

The excitation number of the total system is conserved during the state evolution since the 
excitation number operator N = E^^i(|ei)(ei|+aja^) commutes with the Hamiltonian H. For 
different initial excitation numbers, the system will evolve in different invariant subspaces. 
The zero-excitation subspace is: 

ro = {bii72)|0)i|0)2}, (4) 
and the single-excitation subspace is: 

Ti = {|ei^2>|0)i|0)2,bie2)|0)i|0)2, 

bi^2)|l)i|0)2,bi<?2)|0)i|l)2}. (5) 
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The A^-excitation subspace (A^ = 2, 3, 4, oo) becomes: 

Tat = {^a,^b,^c,^d}, 

Ta = {|ei^2)|iV-l)i|0)2,|ei^2)|iV-2)i|l)2,..., 

|ei(?2)|l)i|A^-2)2,|ei(?2)|0)i|7V-l)2}, 
Tfe = {bie2)|iV-l)i|0)2,|^7ie2)|iV-2)i|l)2,..., 

bie2)|l)i|iV-2)2,|^7ie2)|0)i|iV-l)2}, 
Tc = {bi^2)|A^)i|0)2,bi(72)|A^-l)i|l)2,..., 

M2)|l)l|iV-l)2,|^7l^2)|0)l|iV)2}, 

= {|eie2)|iV-2)i|0)2,|eie2)|7V-3)i|l)2,..., 

|eie2)|l)i|iV - 3)2, |eie2)|0)i|iV - 2)2}. (6) 

Suppose the density operator of the whole system is initially p(0). After an interaction 
time the density operator is given by p(t) = e~*^*p(0)e*-'^*. Taking a partial trace over two 
thermal fields, we can obtain the reduced density matrix for two atoms which is spanned in 
the basis { \e1e2), \e1g2), \gie2)i \9192) }■ Without any initial coherences from both the atoms 
and fields, we calculate the reduced density matrix in each excitation subspace, then sum 
over all the subspaces to get the total reduced density matrix (7). The used Hilbert space is 
cut off at ~ 5(n + 1). For example, when n = 0.1, the Hilbert space is cut off at = 5; 
when n = 1, the Hilbert space is cut off at A^ = 15; when n = 10, the Hilbert space is cut 
oS at N = 60. Therefore, the density matrix Pa has a simple form: 



Pa 



(A OCX 

5^0 

C 

V G* ^ D ) 



(7) 



The values of the elements in Eq. (7) depend on the initial state of the whole system and 
they can be explicitly calculated by numerical simulation. We adopt the Wootters's [34] 
concurrence C{pa) to quantify the entanglement for two atoms, which is defined as 



C(p„)=max|o,yA;-yA^-0;;-yA;}, (8) 

where Ai, A2, A3 and A4 are the eigenvalues arranged in decreasing order of the following 
matrix 

i = Pa{(Ty®ay)pl{ay®ay), (9) 

where p* represents the complex conjugation of pa and ay is the corresponding Pauli matrix. 
Therefore the concurrence for the two atoms takes the form: 

C(pa) = 2max|o, |E| -VAD, |G| - V^C}. (10) 
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It has been shown that using the concurrence to deal with the mixed state of any two qubits 
is an easy evaluation of entanglement of formation, and it satisfies the following important 
conditions of the measure of entanglement: (1) C{pa) = if and only if pa is separable; (2) For 
all local unitary transformations, C{pa) keeps invariant; (3) C{pa) does not increase under 
local general measurements, classical communications, and post-selection of subensemble. 
For unentangled atoms C{pa) = 0, while C{pa) = 1 for the maximally entangled state. 
Although the exact analytic expression for Eq. (10) is very complicated as the excitation 
number becomes infinite, there are several limiting situations where the whole dynamics can 
be simplified. Under these limiting conditions, we can obtain an effective Hamiltonian that 
refiects the essentially physical processes and provides a basic explanation for the numerical 
results. 

3. Generation of Atomic Entanglement 

Assume that the initial state for the two atoms is \eig2)> Then the density matrix of the 
whole system is 

oo 

p(0) = |ei^2)(ei^2| ® J] Pi{ni)P2{n2)\nin2){nin2\. (11) 

Set ni = = ^ = 1, the concurrence is plotted for the cases of large frequency detuning and 
large hopping strength in Fig. 2 and Fig. 3, respectively. In order to explain the underlying 
physics explicitly, we introduce two delocahzed bosonic modes hi and 62, which are defined 
as hi = (ai + a^l\f2 and 62 = (<^i — Gi2)lV^- We first consider the symmetric coupling 
situation, i.e., gi = g2 = g^ and set 5i = 62 = S for simplicity. In terms of the new operators, 
the Hamiltonian H can be rewritten as follows: 

H = 5[h\hi + 5'^hlh2 



hi{S^ + S^) + 62(5+ - ^2^) + H.c. , (12) 



where S[ = {S + J) and S2 = {S — J)^ which indicates that the frequencies of two delocahzed 
field modes are shifted from the cavity frequency due to coherent photon hopping. In the 
interaction picture with respect to Hq = S[h\hi + ^^^2^ we get the atom-field interaction 
Hamiltonian 



e-'^^%i{St + ^2^) + e-'^^%2{St - St) + B.C. . (13) 

Under the conditions 5^ ^ yjfii + \gj\f2 and ^2 ^ \/^2 + ^g/V^i there is no energy exchange 
between the atomic system and the field modes, and the two atoms are coupled to each other 
via exchange of virtual photons. Then the effective Hamiltonian is given by [20,35]: 
^2 2 

Heff = - E + ^4^2)(|e.)(e,| - + A|e,)(e,|] + XiStS:^ + H.c), 

A — 1 1 2 
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(14) 

2 2 / 2 2 / 

where X = ^ + ^ and X = — The evolution of the atomic system is cos(A t)\eiq2) 
— i sin(A't)|5'ie2), which is independent of the field states. We here have discarded the trivial 
common phase factor e~^^ ^ Therefore, for large photon hopping strength the local reso- 
nant thermal fields can induce maximal atomic entanglement, which is distinguished from 
the case for two atoms interacting commonly with one single-mode thermal field [4]. The 
eflFective Hamiltonian (14) is valid under the condition that the difference between 5 and J 
is much larger than g. When 5 ^ the eigenfrequency of the normal delocalized mode 62 
approaches the atomic transition frequency so that the atoms can exchange energy with this 
delocalized mode and the eflFective Hamiltonian (14) is invalid. This leads to the atom-cavity 
entanglement and deteriorates the atom- atom entanglement, which accounts for the appear- 
ance of the dips in Fig. 2 and Fig. 3 for 6 near J. In the regime 5 ^ J ^ the atom-cavity 
system undergoes Rabi oscillations, and each atom is not only entangled with each other 
but also with the cavity modes. One can maximize the atom-atom entanglement by choosing 
the optimal ratio of g to and t. The entanglement is very sensitive to the ratio S^/ g- This 
accounts for the maxima and abrupt changes of the concurrence inside the dips of Fig. 2 and 
Fig. 3. Due to the spread of the Rabi frequencies corresponding to diflFerent photon num- 
bers, the atoms cannot be completely disentangled with the cavity modes and the maximal 
atom- atom entanglement cannot be obtained. In fact, the physics discussed here is similar 
to that discussed in Ref. [35] which proposed the generation of maximally entangled states 
for two atoms via dispersive interaction with a single cavity mode. 

When two atoms are both initially in their excited states |eie2), we find that the atomic 
entanglement could also be induced by two thermal fields with appropriate parameters, 
as shown in Fig. 4. In the case of Fig. 4(a), the atomic entanglement arises from either 
atom's emission of one photon into the same cavity modes. When the probability for this 
event is zero, the entanglement vanishes. In Fig. 5, we plot the expectation value of total 
atomic excitation number {a^i + <^z2) versus gt for the same parameters as Fig. 4(a). The 
result shows the peak of concurrence coincides with the case when (a^i + (7^2) is nearest to 
zero, which corresponds to the regime where probability for either atom's emission of one 
photon is maximized. Interestingly, for symmetric atom-cavity couplings the entanglement 
exhibits periodical sudden birth and death when the atomic transition is slightly detuned 
from at least one delocalized mode, as shown in Fig. 4(a). Note that the atoms start from 
a product state, which is distinguished from the previous studies of entanglement sudden 
birth and death in which the atoms start from an entangled state [28-30]. The duration of 
the entanglement death is much longer than that of the entanglement life. Unlike previous 
examples, the peaks of the concurrence do not decrease monotonically as the interaction 
time increases. The entanglement dynamics for the initial state \g1g2) is similar with that 
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for initial state |eie2), as shown in Fig. 4(b). The entanglement dynamics for both the initial 
states 16162) and \g1g2) is similar to the case when the two atoms interact with a common 
cavity mode [4,5]. On the other hand, for asymmetric atom-cavity couplings, i.e., gi ^ 5^2, 
the atomic entanglement can be generated when the atoms resonantly interact with two 
thermal fields, as shown in Fig. 4(c). The vanishes of entanglement in Fig. 4 corresponds to 
entanglement sudden death that also appears in the system with two atoms independently 
interact with the respective cavities [28] . The considered system involves four subsystems so 
that the whole system evolution is very complex and does not exhibit periodical features. 
This accounts for the fact that distributions of peaks and dips in Fig. 3 and Fig. 4 are 
somewhat random. 

4. Freeze of Atomic Entanglement 

Moving forward, we also investigate the system evolution when the atoms are initially in the 
maximally entangled state \^rn) — {\^ig2) + \gi^2)) / V^- In the resonant interaction situation, 
as shown in Fig. 5, for J 5^, the entanglement can be frozen. This is due to the fact in 
this case the conditions 5'-^ ^ \/n\ + \gj\f2 and ^ + ^glV^ are satisfied, so that 
the effective Hamiltonian of Eq. (14) dominates the system evolution and the maximally 
entangled state \^rn) becomes an eigenstate of this effective Hamiltonian. As the values of 
J/g decreases, the concurrence decreases. Especially when J/g ^ 0^ the sudden death of 
entanglement emerges, since in this case each atom independently interacts with a thermal 
field, which causes the leakage of the atomic coherence into the fields. 

The infiuence of temperature, i.e., mean photon numbers n, on the atom- atom entangle- 
ment are taken into consideration, as plotted in Fig. 7. The result shows that the larger the 
mean photon number, the smaller the concurrence. This phenomenon can be understood 
in the following way. The probability that the atoms exchange energy with the fields in- 
creases with mean photon numbers. When the photon numbers are large, the conditions S[ 
^ \/n\ + \gj\f2 and ^ + ^g/V^ are not satisfied and the effective Hamiltonian 
(14) is not valid. When the photon number is sufficiently small so that the effective Hamil- 
tonian (14) is valid, the effective atom-atom coupling strength is A' = g'^ / J = O.lg^ which 
is photon- number-independent. When X't = 7r/4, i.e., gt = 2.57r, the excitation is approxi- 
mately equally shared by the two atoms and the concurrence reaches the maxima. At the 
time t = 7r/(2A') = 57r, the excitation is completely transferred from the first atom to the 
second one. So that the concurrence vanishes. 

5. Atomic Entanglement Dynamics with Dissipation Being Included 

In all the above discussions, we have just supposed that the entire system is ideally isolated 
from the outside environment. Taking cavity decay and atomic spontaneous emission into 
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account, the master equation for the density matrix p{t) of the system can be expressed as: 
Pit) = -t[H,pit)] 

2 

+ -(n + 1) ^ [2aip{t)al - a\aip{t) - p{t)a\ai] 
1=1 

2 

^ i=i 

2 

+|(^ + 1) E {^S-p{t)St - St S- pit) - pit)sts-] 

1=1 

2 

+lnY,[^Stp{t)Sr - SrStp{t) - p{t)SrSti (15) 
^ i=i 

where k and 7 denote the cavity decay rate and the atomic spontaneous emission rate, re- 
spectively. The dependence of the concurrence on the cooperative parameter (7^^ = 5'^/(/^7) 
is plotted, as shown in Fig. 8. When the detunings between atomic transition and the de- 
localized modes are much larger than the atom-cavity coupling strength, the entanglement 
dynamics is much more sensitive to the atomic decay than to the cavity loss as in this case 
the cavity modes are only virtually excited and they are decoupled from the atomic entan- 
glement dynamics. Therefore, the scheme is favorable when 7 <C ^. Here we choose 7 = 0.1/^ 
in Fig 8(a) and (c), and 7 = in Fig 8(b) and (d). The difference between Fig 8(a) ((c)) and 
Fig 8(b) ((d)) is that the concurrence decreases quickly as the atomic decay increases, which 
demonstrates the cavity decay is not the dominant factor deteriorating the entanglement. 

6. Conclusions 

In conclusion, we have studied the entanglement dynamics of two identical atoms interact- 
ing with two single-mode thermal fields based on the coupled-cavity system. The atom-atom 
concurrence is calculated under various conditions. The results show that the entanglement 
behaviors depend upon the photon hopping strength, the atom-cavity detuning, and the 
initial atomic state. It is demonstrated that under the situation where each atom is resonant 
with the local field mode, the maximal entanglement for two atoms can be generated with 
one atom in the ground state and the other in the excited state initially, which is impossible 
for the case with two atoms interacting with a single thermal cavity mode. When the two 
atoms are both initially in the excited state, the entanglement sudden birth and death oc- 
cur periodically, and the peaks of the concurrence do not decrease monotonically when the 
evolution time increases. Furthermore, we show that the initial two-atom maximal entangle- 
ment can be frozen with suitable choice of the parameters. The eflFect of dissipation on the 
entanglement is analyzed. 
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List of Figure Captions 

Fig. 1. (Color online) The proposed experimental setup. Two two-level atoms are respectively 

trapped in two coupled cavities, each of which is initially in a thermal state. 

Fig. 2. (Color online) Atom- atom concurrence as a function of the evolution time and photon 

hopping strength with n = 1 and S = lOg. The atoms are initially in the state \e1g2). 

Fig. 3. (Color online) Atom-atom concurrence as a function of the evolution time and photon 

hopping strength with n = 1 and J = 25g. The atoms are initially in the state \e1g2). 

Fig. 4. (Color online) Atom-atom concurrence as a function of the evolution time and hopping 

strength when n = 0.1 : (di) gi = g2 = g for the initial state \e1e2)] (b) gi = g2 = 9 for the 

initial state \g1g2)] (c) S = 0^ gi = J and g2 = g for the initial state \e1e2). 

Fig. 5. (Color online) Dash red line represents the atomic population inversion (a^i + (7^2) 

and solid black line represents the average photon number (a|ai), which are as a function of 

the evolution time when n = 0.1, J = 205^ and S = 18.5g for the initial state \e1e2). 

Fig. 6. (Color online) Atom-atom concurrence as a function of evolution time and hopping 

strength with S = and n = 0.1. The atoms are initially in the entangled state {\e1g2) + 

\9ie2))/V2. 

Fig. 7. (Color online) Atom-atom concurrence as a function of evolution time with 5 = 0, J 
= lOg. The atoms are initially in: (a) \e1g2)] (b) {\e1g2) + \gie2))/V^- 

Fig. 8. (Color online) Atom-atom concurrence as a function of the evolution time gt and 
cooperative parameter (7^^ when n = 0.1. The atoms are in the initial state \e1g2) with 
6 = and J = lOg: (a) 7 = 0.1^; (b) 7 = While the atoms are in the initial state 
(1^192) + 191^2))/ with S = 15g and J = 5g: (c) 7 = 0.1/^; (d) 7 = ^. 
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